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Abstract. We show that the only orthogonal polynomials with a generating function of the form 
F (xz — az 2 ) are the ultraspherical, Hermite, and Chebyshev polynomials of the first kind. The 
generating function for the Chebyshev case is non-standard, although it is easily derived from the 
usual one. 



1. The question 

Hermite polynomials H n {x) are one of the most important families of orthogonal polynomials in 
mathematics, appearing in probability theory, mathematical physics, differential equations, combi- 
natorics, etc. One of the simplest ways to construct them is through their generating function, 



—H n {x)z n = exp (xz - z 2 /2) 



n=0 



On the other hand, Chebyshev polynomials of the second kind U n (x) are another important fam- 
ily of orthogonal polynomials (appearing in numerical analysis, for example), with a generating 
function 



Y,Un{x)z n = - l — -. 

^— ' 1 — XZ + z z 

n=0 



Note that both of these functions have the form F (xz — az ), with F(z) = e z , respectively, F(z) = 
1 

1-2 - 

Question 1. What are all the orthogonal polynomials with generating functions of the form 

F (xz — az 2 ) 

for some number a and function (or, more precisely, formal power series) F? 

A reader interested in further context and background for this question may want to start by reading 
Section 



2. The method 

The following is the first fundamental theorem about orthogonal polynomials. 

Theorem (Favard's theorem). Let {Pq(x), Pi(x), P2(x), . . .} be a monic orthogonal polynomial 
family. That is, 
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Each P n is a polynomial of degree exactly n (polynomial family). 
Each 

Pn{x) = x n + lower order terms 

(monic). 

For some non-decreasing function M with lim a ._ > _ 0O M(x) = 0, Hindoo M (x) = 1, and 
infinitely many points of increase ( equivalently, for some probability measure supported on 
infinitely many points) the Stieltjes integral 



P n {x)P k {x) dM(x) = 



for alln k (orthogonal). 



Then there exist real numbers {(3 Q , (3 2 , ■ ■ ■} and positive real numbers ^2, <^3, • • •} such that 
the polynomials satisfy a three-term recursion relation 

xP n = Pn+l + fi n Pn + ^nPn~l 

(to make the formula work for n = 0, take P-\ = 0). 



Actually, Favard's theorem also asserts that the converse to the statement above is true: if we 
have a monic polynomial family satisfying such a recursion, these polynomials are automatically 
orthogonal for some M(x), and the case of M with only finitely many points of increase can also 
be included. See HChi78l, Chapter 1] or HIsm051 Chapter 2] for the proof. We only need the "easy" 
direction stated above, which we now prove. 



Proof. Every polynomial P of degree A; is a linear combination of {Pi : < i < k}. Therefore P 
is orthogonal to all P n with k < n. Since 

xP n (x) = x n+1 + lower order terms, 

we can expand 

xP n {x) = P n +l{x) + C n)n P n (x) + C ni „_iP„_i (x) + . . . + Cn^P^x) + C nfi P (x) 

for some coefficients c ni „, . . . , c„ j0 . On the other hand, for k < n — 1 orthogonality implies 



P k (x)P k (x)dM(x) = j (xP n (x))p k (x)dM(x) = j P n (x)(xP k (x))dM(x) = 



since deg (xP k ) < n. Since M has infinitely many points of increase, P k cannot be zero at all of 
those points, and as a result 



oo 

2/ 



P k \x)dM(x) > 0. 
Therefore c n)k — for k < n — 1, so in fact 

xPn^X*) -P n _|_i(x) -|- (' n ji P n yX) -\- C nn —\P n —\(jX^) 
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only. Denoting (3 n = c n n and u n = c^ n _i, we get the formula. It remains to note that 

/*oo i*co . 

C-. n _l / P n - 1 (x)P n -l(x)dM(x) = / P n - 1 (x)(xP n (x) - P n +\(x) - Cn.nPJx) | il,\t(.r) 



=n,n-i / P„-i(a;)P n _i(a;) dM(x) = j P n - X {x) (xP n (x) - P n+1 (x) - c n ,„P n (x)) dM(c 

J — CO J —CO 

P n -x(x)(xP n (x)^dM(x) 
/ xP n . l {x))p n {x)dM{x) 



—oo 
oo 



P n (x)P n (x)dM(x). 

Since both P^_ 1 (x) dM(x) and P^{x) dM(x) are positive, so is u n = c ni7l _i. □ 

3. Examples 

Good references for polynomial families are Wikipedia and [KS98|. Beware that in the following 
examples, we use for monic polynomials the notation which usually appears for other normaliza- 
tions, so our formulas may differ from the references by a re-scaling. 

Example 1. The orthogonality relation for the Hermite polynomials is 

(1) H n {x)H k {x) (-^=e- x2 ' 2 \ dx = 0. 



Thus they are orthogonal with respect to the normal (Gaussian) distribution. They satisfy a recursion 

xH n (x) = H n+ i(x) + nH n _ x {x). 
Example 2. The orthogonality relation for the Charlier polynomials is 

(2) Y,C n {i)C k {i)(e- x -\ =0, 

so that 

W 1 

i=0 

which is sometimes also written as 

dM(x) = e" 1 V -5 l (x). 
L — ' %\ 

i=0 

Thus the Charlier polynomials are orthogonal with respect to the Poisson distribution. They satisfy 
a recursion 

xC n [x) = C n+ i(x) + (n + 1) C n (x) + n C n -i(x). 

Example 3. The Legendre polynomials are the family of orthogonal polynomials a student is most 
likely to encounter in an undergraduate course. In a linear algebra course, one sees them in the 
applications of the Gram-Schmidt formula, since their orthogonality relation is simply 

P n (x)P k (x)dx = 0. 

l 
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In a differential equations course, one sees them as a solution of the Legendre equation. This 
equation, in turn, arises after the separation of variables in the heat equation or Laplace's equation 
in three dimensions, or on a sphere. The recursion satisfied by the Legendre polynomials is 

xP n (x) = P n+ x{x) + — — -P n _i(x) 
An z — 1 

Example 4. The ultraspherical (also called Gegenbauer) polynomials are orthogonal with 
respect to 

M(x) = (1 -x 2 ) A ^ 

(with integration restricted to [—1, 1]) for A > — |, A 0. Note that the Legendre polynomials are a 
particular case corresponding to A = \. Just as the Legendre polynomials are related to the sphere, 
the other ultraspherical families for half-integer A are related to higher-dimensional spheres. They 
satisfy a recursion 

(3) xC^\x) = C^l{x) + A{n [^_ 1){n l x) C ( n\(x), 

see Section 1.8.1 of [KS98|. Another important special case are the Chebyshev polynomials of the 
second kind U n , which correspond to A = 1 and are orthogonal with respect to y/l — x 2 . 

For A = 0, the coefficient u\ as written in formula Q is undefined. The polynomials orthogonal 
with respect to ^=f are the Chebyshev polynomials of the first kind T n , for which the recursion 
@ holds for n > 2 (with A = 0), but 

xT x (x) = T 2 {x) + ^T (x). 

Remark 1. Suppose {P n } form a monic orthogonal polynomial family, with orthogonality given 
by a function M and the recursion 

(4) xP n = P n+l + /3 n P n + w n P n _i. 
Then for r ^ 0, the polynomials 

Q n {x) = r n P n (x/r) 
are also monic, orthogonal with respect to N(x) = M(x/r), and satisfy 

xQ n = Qn+i + r(3 n Q n + r 2 uj n Q n ^ X - 

Indeed, from equation Q, 

r n+ \x/r)P n {x/r) = r n+1 P n+1 (x/r) + f3 n r n+1 P n (x / r) + u n r n+1 P n ^(x/r), 
which implies the recursion for {Q n }. 

4. The result 

Theorem 1. Let a > and F(z) = J2^Lo c nZ n be a formal power series with cq — 1, c\ — c ^ 0. 
Define the polynomials {P n : n > 0} via 

oo 

(5) F (xz - az 2 ) = c n P n (x)z n 

71=0 

(ifcn = 0, P n is undefined). These polynomials form an orthogonal polynomial family (which is 
automatically monic) if and only if 



A CHARACTERIZATION OF ULTRASPHERICAL POLYNOMIALS 

{P n } are re-scaled ultraspherical polynomials, 

for X > — \, X 7^ 0, and b > 0. In this case 

F(z) = 1 + ' / ' 




4a X 



Xb \(l-bz) x 

for c 7^ 0. The choice c = Xb gives simply F(z) = ( . 1 _^ z - )A . 
{P n } are re-scaled Chebyshev polynomials of the first kind, 



P n (x) = T n I x 




for b > 0. In this case 



4a 



F(z) = 1 + \ In 1 ' 



b W-bz 



for 0. The choice c = b gives simply F(z) = 1 + In (jz^)- 
{P n } are re-scaled Hermite polynomials, 




for a > 0. In this case 



P n (x) = H„ ( \ x 
la 



F(z) = 1 + - (e az - 1) 
a 



for 0. The choice c = a gives simply F(z) = e° 
Remark 2. If equation Q holds, so that 



1 + c n (xz — az 2 ) n = F [xz — az 2 ) = 1 + c n P n (x)z n , 

n=l n=l 

then clearly also 

oo oo 

1 + Yl CCn ( XZ - = F ° ( XZ - ^ = 1 + Yl Cc nPn{x) Z n 

n=l n=l 

for any C ^ and Fc(z) = 1 + C(F(z) — 1). This is the source of the free parameter c in the 
theorem. 
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Proof of the Theorem. Using the binomial formula, expand 



F (xz — az 2 ) = c n (xz — az 2 ) n 

n=0 

oo oo 

= c n x n z n - c n n[xz) n ~ x az 2 + . . . 

n=0 n=l 

oo 

= 1 + ^{c n x n - a{n - l)c n _ 1 x n - 2 ^z n + ... 



n=l 



If for some n > 2, c n = while c n _i 7^ 0, then comparing the coefficients of z n in the preceding 
equation and expansion (|5]), we see that the coefficient is non-zero on the left and zero on the 
right. So all c n 7^ 0, and we may denote d n = Then the same coefficient comparison gives 
Pq(x) = 1, and forn > 1 

p n ( x ) = x n - a{n - l)d~ l x n ~ 2 + ... 

Using this equation for both n and n + 1, we then get 

x p n = _ a (n - lX 1 ^" 1 + . . . 

= P n+ i - a(n - l)^ 1 ^"" 1 + and~\ l x n ~ 1 + ... 

If we want the polynomials to be orthogonal, by Favard's theorem they have to satisfy a three-term 
recursion relation 

xP n — Pn+l + finPn + 

(note that {P n } are clearly monic). We see that (3 n = 0, and 

ui n = a(nd~\ x - (n - l)d~ x j 

for n > 1 . 

Now expanding further in the binomial formula, 



OG 



c n (xz - az 2 ) n = c n x n z n - c n n(xz) n 1 az 2 + y^c n -^ — -{xz) n 2 (az 2 ) +... 

n=0 n=0 n=l n=2 

= 1 + C\xz + (c 2 x 2 — aci) z 2 + (c 3 x 3 — 2ac 2 x) z 3 

+ Ux n - a(n - l)c n . lX n - 2 + a 2 {n ~ 2) ^~^ c n . 2 x n A z n + . . 

Thus 

Po(x) = 1, P\{x) = x, P 2 (x) = x 2 — ad 2 x , Pz(x) = x 3 — 2ad 3 l x 
and for n > 4, 

p n ( x ) = x ™ _ a (n - l^V" 2 + a a ( "" 2)( "" 3) rfXV' 4 + • • • 
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Therefore for n > 4, 

/ 2 (n-2)(n-3) x x 2 (n - \){n - 2) x x 
xP n — P n +i — u n P n -i — I a d n d n _ 1 — a d n+l d n 



- a((n - IX 1 - nd-l^ain - 2)Ci)^ 3 + 



a formula which also holds for n = 3. For this to be zero we need 



(71 — 1) / \ 

d n+ i — d n -i - \ {n - l)d n+ i - nd n J = 0. 



2 

or 

(n-3) J (n-1) 



Thus for n > 3, 

(n + l)d n+ i = 2nd n - (n - l)d n _i. 
The general solution of this recursion is 

nd n = a + b{n — 1) 

for n > 2. Since all d n ^ 0, a,b cannot both be zero. Therefore 

r _ a +6(n-i) nffi^L _ nr=i(«+ao „ 

n n\ n\ 

for n > 2 and 

n=2 

If a ^ 0, 6 ^ 0, then 



(6) 



n=l 



1 + - ((1 - bz)- a / b -!)=! + - f 1 - 1 
a vv 7 7 a \(1- fe) a / b 



If a = 0, b ^ 0, then 

00 Wl— l / 1 \ 

(7) F W = H^ T2 -1 - -°h.<l-ta) = 1+ 5l„ _ , 

n=l N ' 

which can also be obtained from the preceding formula by using L'Hopital's rule. Finally, if a 
6 = 0, then 

00 n-l 

(8) F(z) = 1 + z n = 1 + -(e az - 1). 

n=l 

Moreover, 

(n- 1)6 + 2a 



oo n = an-. 



((ra- l)6 + a)(n6 + a)' 
Since for orthogonality, we need w n > 0, clearly 6 > 0. If 6 = 0, then 

2a 

u n = — n > 
a 
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as long as a > 0. The polynomials with this recursion re- scaled Hermite polynomials. We recall 
[KS98, Section 1.13] that the generating function for standard (monic) Hermite polynomials is 



—H n {x)z n = exp (xz - z 2 /2) 



n=0 



which is of the form ([5]) with F(z) = e z , and the generating function ([8]) is obtained from it by a 
re-scaling and a shift from Remark [2] 

If b > 0, a 7^ 0, we denote A = a/b and get 

a n(n + 2A - 1) 



OJr. 



Since 



we have A > — 1. Since 



b {n + A - l)(n + A) ' 
a 2 



uj 2 



fel + A' 
a 2(1 + 2A) 



b (1 + A)(2 + A)' 

we have moreover A > —\. It is now easy to see that this condition suffices for the positivity of 
all ui n ; indeed, the corresponding polynomials are re- scaled ultraspherical polynomials. We recall 
[KS98, Section 1.8.1] that the generating function for standard (monic) ultraspherical polynomials 

2«nr=o(A-OwA)^n_ i 



n\ " * )Z ~ (l-2xz + z^ 

which is of the form (|5]) with F(z) = ji_\z) x ' anc ^ me g enerat ing function ([6]) is obtained from it by 
a re- scaling and a shift from Remark [2j 

Finally, if b > 0, A = a = 0, then 

a 

UJ n = ~r 

b 

for n > 2, but u>i = 2^. These are precisely recursion coefficients for the re-scaled Chebyshev 
polynomials of the fist kind. The standard generating function [KS98, Section 1.8.2] for (monic) 
Chebyshev polynomials of the first kind is 

OO ^ 



1 — 2xz + z 2 ' 

n=0 

so it is not of the form ([5]). However, 

n — 1 V * 



n=l 

Term-by-term integration with respect to z gives 



1 — 2xz + z 2 



C + V — TJx)z n = -- In (1 - 2xz + z 1 ) = - In ( — 1 — ^ . 

n=l x v it/ 

with C = 0, which is of the form (|5]) with F{z) = | In (tz^-). The generating function @ is 
obtained from it by a re- scaling and a shift from Remark [2} □ 
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5. The history 



The question of characterizing various classes of orthogonal polynomials has a long and distin- 
guished history, see HAS90H for an excellent survey up to 1990. The study of general polynomial 
families goes back to Paul Appell in 1880 [App80|, who looked at polynomials with generating 



P n {x)z n = A{z) exp(xz) 



functions of the form 

oo 1 

Y- 

71=0 

for some function A (z). These are now called Appell polynomials. Later, they were generalized to 
Sheffer families with generating functions 



(9) 



oo 1 

V -Mx)z n = A(z) exp(xU(z)) 



n=0 



for some functions A(z), U (z). The prototypical "orthogonal polynomials characterization result" 
is Meixner's 1934 description of all orthogonal polynomials with the Sheffer-type generating func- 
tions [M ei34L On the other hand, among Appell polynomials, only Hermite polynomials are or- 
thogonal. See Figure 1 . 




Figure 1 . The relationship between the class of orthogonal polynomials and Shef- 
fer, Appell and Meixner classes. "H" stands for Hermite. 



Example 5. One way to state Meixner's result is that orthogonal polynomials with generating 
functions Q satisfy a three-term recursion 

xP n (x) = P n+ x{x) + {na + /3 ) P n (x) + n{{n - l)b + 1) P n - X {x), 

for some a, b, (3 (up to re-scaling). The function M a b / g for which these polynomials are orthogonal 
can be written down explicitly, but for different values of the parameters these functions look quite 
different. For example, for a = b = (3q = 0, we get the Hermite polynomials, with a continuous 
orthogonality relation ([T]). On the other hand, for a — /3q — 1, b — 0, we get the Charlier poly- 
nomials, with a discrete orthogonality relation (|2]). Other polynomials in the Meixner class carry 
the names of Laguerre, Krawtchouk, Meixner, and Pollaczek, and are orthogonal with respect to 
gamma, binomial, negative binomial, and Meixner distributions. 



Besides nice generating functions, the Meixner class has many other characterizations and applica- 
tions, see [DKSC08] for an excellent (but advanced) survey. Perhaps for this reason, many general- 
izations of this class have been attempted. The most popular of these are probably the g-deformed 
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families. One approach (there are several) extends the Sheffer class by looking at the generating 
functions of the form 



(1 -q)U{q k z)z 

(after appropriate normalization, one gets the Sheffer form for q — > 1). In this case the analog of the 
Meixner class are the Al-Salam and Chihara polynomials [ASC87]. For the study of two different 
types of g-Appell polynomials, see HAS67[|AS95L 

A different generalization of the Sheffer class are generating functions of the general Boas-Buck 
HBB64II type: 

oo 

Y,c n P n (x)z n = A(z)F(xU(z)) 

n=0 

for F(z) = Xl^Lo c «- 2 ™ w ^ tn Co — 1- The usual case corresponds to F(z) = e z . In the Boas-Buck 
setting, the problem of describing all orthogonal polynomials is wide open. The Appell-type class 
(with U (z) = z) in this case consists of the Brenke polynomials, and at least in that case all the 
orthogonal polynomials are known [Chi68]|. 

Now note that in the Sheffer/Meixner case in equation (|9]), corresponding to F(z) = e z , the gener- 
ating function has an alternative form 

—PJx)z n = A(z)exp(xU(z)) = exp(xU(z) + log A(z)). 
/ — ' n! 

n=0 

So another interesting class to look at are (all or just orthogonal) polynomials with generating 
functions 

oo 

Y,CnPn{x)z n = F{XU{Z) - R{z)), 
n=0 

which again gives the Sheffer/Meixner families for F(z) = e z . 

The case F(z) = j^- appears in Free Probability [NS06J, see Section 3 of HAns03H for the author's 
description of the "free Meixner class", which is in a precise bijection with the Meixner class 
(except for the binomial case [BB06]). Here again, one can write the generating function in two 
ways: 

1 1 

A(z) — = 7 

K J l-xU(z) i (~v(*) 

1 \ X A(z) A(z) 

More generally, Boas and Buck proved the following result. 
Theorem. [BB561 The only functions F with F(0) = 1 such that 
(10) A(z)F(xU(z)) = F(xu(z) - R(z)) 

are F(z) = e z and F(z) = n^jx for some A. 

So as an alternative to the Boas-Buck formulation, we are interested in orthogonal polynomials 
with generating functions of the form F(xu(z) — R(z)), or at least in the Appell-type subclass 
F(xz — R(z)). For general R, even this seems to be a hard question. However, the orthogonal 
Appell polynomials are only the Hermite polynomials, with the exponential generating function 

exp [xz — z 2 /2j . 
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On the other hand, the orthogonal free Appell polynomials are only the Chebyshev polynomials of 
the second kind, with the ordinary generating function 

1 

1 — (xz — z 2 ) 

Moreover, R(z) = az 2 appears naturally in combinatorial proofs of the usual, free, and other 
central-limit-type theorems (see for example Lecture 8 of [NS06|). Thus it is reasonable to consider 
F (xz — az 2 ) first, which leads to Question [I] Conversely, the answer to that question indicates 
that interesting generating functions (and also, potentially, interesting non-commutative probability 
theories) arise precisely for F covered by the Boas-Buck theorem above, plus in the exceptional 
case F(z) = 1 + log y^- not covered by that theorem. On the other hand, see HASV861 lDem09B for 
some negative results. 

Acknowledgments. I am grateful to Harold Boas for a very careful reading of the paper, and for 
many valuable comments, including Remark[2j 
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